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Q ' ABSTRACT 
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We report on the discovery of strongly variable Fe II and Mg II absorption 
lines seen at z — 1.48 in the spectra of the z = 4.05 gamma-ray burst (GRB) 
q ■ 060206 obtained between 4.13 to 7.63 hours (observer frame) after the burst. In 

particular, the Fe II line equivalent width (EW) decayed rapidly from 1.72 ± 
0.25 A to 0.28 ± 0.21 A, only to increase to 0.96 ± 0.21 A in a later spectrum. 
The Mg II doublet shows even more complicated evolution: the weaker line of the 



> 



doublet drops from 2.05 ± 0.25 A to 0.92 ± 0.32 A, but then more than doubles 
^ ! to 2.47 ± 0.41 A in later data. The ratio of the EWs for the Mg II doublet is also 



variable, being closer to 1 : 1 (saturated regime) when the lines are stronger and 
becoming closer to 2 : 1 (unsaturated regime) when the lines are weaker, consistent 
with expectations based on atomic physics. We have investigated and rejected 
the possibility of any instrumental or atmospheric effects causing the observed 
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strong variations. Our discovery of clearly variable intervening Fe II and Mg II 
lines immediately indicates that the characteristic size of intervening patches of 
Mg II "clouds" is comparable to the GRB beam size, i.e., about 10 16 cm. We 
discuss various implications of this discovery, including the nature of the Mg II 
absorbers, the physics of GRBs, and measurements of chemical abundances from 
GRB and quasar absorption lines. 

Subject headings: galaxies: ISM, gamma rays: bursts, intergalactic medium 



1. Introduction 



The optical spectrum of a gamma-ray burst (GRB) afterglow is generally a power-law 
with absorption lines superimposed. These absorption lines correspond to intervening mate- 
rial between the source and the observer. Typically, the highest-redshift li nes are thought to 
be associated with the host galaxy (e.g., see discussion in lJha et al. 11200 ll ). and thus provide 
a measurement of the redshift of the burst. There might be several lower- redshift absorption 
lines caused by intervening systems. Spectra of luminous GRB afterglows have opened a new 
win dow on absorption-lin e studies along lines of sight that are not associated with quasars 
(cf. IVreeswijk et al.ll2003l ). 



Recent studies of GRB absorption spectra have produced an unexpected result. iProchter et al. 



( 120061 ) reported statistically significant evidence for a much higher incidence of strong Mg II 
absorbers along GRB sight lines than that expected from studies of quasars of comparable 
redshift. They considered and rejected three possible scenarios to produce such an effect, 
namely dust associated with the Mg II absorbers, ejection of Mg II by the GRB, and strong 
gravitational lensing by galaxies associated with the emission. 



Frank et al. I (120061 ) proposed a simple geometric solution to the quandary, with resulting 



predictions that are readily testable. Namely, if the Mg II absorbers seen in both GRB and 
quasar spectra are truly intervening, then the Prochter et al. result could be caused by 
different sizes of the GRB and quasar emission regions. If the Mg II absorbers have physical 
sizes comparable to the emitting regions of quasars and GRBs, then a background source 
with a larger emission region would "see" a smaller effective absorption column density than 
a background source with a smaller emission region (see Frank et al. for details of their 
proposed scenario). 

One immediate prediction of the Frank et al. scenario is that for strongly dynamic 
objects like GRB afterglows, the strength and structure of intervening absorption lines in 
GRB spectra should vary in some, if not all, cases. Motivated by this prediction, we have 
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examined our existing, multi-epoch spectroscopic data of GRB 060206 and have discovered 
such a behavior. In this paper, we provide the first observational evidence of this prediction. 
In fact, as far as we know this is the first report ever of clearly evolving intervening (i.e., non- 
intrinsic) ground state absorption lines seen in spectra of any cosmological object. There is 



evidence for variation in fine structure lines in the GRB host environment (cf. lVreeswijk et al. 



20061 ; iDessauges-Zavadsky et alll2006l ). 



GRB 060206 triggered Swift-B A T on 2006 Feb 6 at 4:46:53 UT (IMorris et al. 1120061 ). 
A likely afterglow was identified by iFynbo et al. I (120061 ). who also determined that the 
afterglow was at high redshift (z = 4.048 ). The evolution of a bright and complex after- 
glow was reported by sev eral groups, and I Wozniak et al.l (120061 ). IStanek et al. I (120061 ). and 
Monfardini et al. I (120061 ) have published detailed photometric data for this "anomalous" 



event. The burst also showed complex behavior in the X-rays as seen with Swift-XBT in- 



strum ent, but the overall X-ray evolution follows that seen in the optical (IStanek et al. 
2006h . 



2. Spectroscopic Observations 



The brightness of the optical transient (OT) associated with GRB 060206 allowed us to 
take spectra for seve ral hours after the bu rst. Spectra of the OT were obtained with the 
FAST spectrograph (IFabricant et al. II1998T ) at the Cassegrain focus of the 1.5m Tillinghast 
telescope at the F. L. Whipple Observatory at 08:25, 09:02, 09:33, 10:05, 10:41, 11:37, 
12:09 (all times UT). The OT was visible on the telescope acquisition camera and could 
be identified and placed on the spectr ograph slit for the observations. Each spectrum was 
taken at or near the parallactic angle (jFilippenko Ill982l ). The spect ra were reduce d in the 
standard manner with IRAF0 using the optimal extraction package of iHorne I (I1986I ). Night- 
sky emission lines were removed by fitting a line to two sky regions on either side of the OT 
column-by-column; the value determined from the linear function at the position of the OT 
extraction was then subtracted. Wavelength calibration was provided with an HeNeAr lamp 
taken immediately after each OT spectrum. Minor adjustments to the wavelength solution 
were made based upon night-sky lines in the OT frames. There were clouds during some 
of the exposure, but the signal-to- noise ratio (S/N) is comparable for all the spectra. The 
resolution of each spectrum is 6.2 A as measured from the FWHM of night-sky emission 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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lines. 



We flux-calibrated the spectra with spectrophotometric standard stars using our own 
software. Relative fluxes in our OT spe ctrum are accurate to 4% over the observed wave - 
length range. Following the arguments of lWade fe Hornd (119881 ) and lMatheson et al. I (120001 ). 
telluric correction was accomplished by dividing the GRB spectra with a normalized spec- 
trum of the flux standard, which was set to unity for regions unaffected by atmospheric 
features. For our spectra, telluric correction was only applied over the regions 6250-6360 A 
and 6840-7410 A. We scaled the strength of the absorption b y the difference of the airmasses 
to the 0.6 power to correct the saturated B band properly (IWade fe Hornd Il988l ). Most of 
the correction is 2% or less, with a 5% correction in the 6250-6360 A region and 25% in 
the strong region of the B band. 

In the S/N weighted average spectrum, shown in Figure (U a damped Lyman-alpha 
absorption system (DLA) is present at A = 6145 A. A number of strong and weak absorption 
metal line systems are su perimposed on the continuum, identified as systems at z = 1.48 
and z = 4.05 (reported by iFynbo et al. II2006I ). The very strong DLA system at A = 6145 A 
indicates that the other absorbers at z = 4.05 most likely represent the host galaxy of the 
GRB. In the intervening metal line system at z — 1.48, our data show absorption from the 
transitions Fe II A2600 and Mg II AA2 796, 2803. These are fairly typical in absorption line 



systems (e.g., ISteidel fe Sargent Ill992l ). Two absorption lines at A = 6331.3 A and 6338.3 A 



represent a Si IV system at z 
5525.9 A. 



3.54, with a strong corresponding Lya absorption line at 



3. Spectroscopic Monitoring 

Motivated by the prediction of Frank et al. (2006) that the absorption lines along GRB 
sight lines should vary in some cases, we examined our time series of GRB 060206 spectra 
to look for this effect. Variation of the intervening z = 1.48 Fe II and Mg II lines is obvious 
even upon casual inspection of the spectra, as shown in Figure |5J These lines are seen to 
vary in strength, profile, and possibly velocity. We have checked the original 2-D spectra 
and no cosmic ray features are apparent that would cause this variation. In addition, night- 
sky emission lines and telluric features were carefully removed with several variations in 
technique to ensure that reduction processes did not introduce the observed variability. In 
fact, there are no sky lines or telluric features at all near the strongly variable z = 1.48 Fe II. 

To measure the magnitude of variations in the equivalent widths (EWs) of the z = 1.48 
absorption lines, we fit each of the Mg II lines with multiple Gaussians and the Fe II line with 
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a single Gaussian. The error of the EWs are calculated from the conservative fit uncertainties. 
The resulting EW measurements are listed in Tableland plotted in Figure [3) The first three 
spectra show a gradual decrease in the EW of the Mg II A2796 line accompanied by a sharp 
decline in the EW of the weaker Mg II A2803 line. As a result, the ratio of the doublet 
increases from 1.2 : 1 in the first spectrum to about 2 : 1 in the third. This is followed in 
the fourth spectrum by a sharp rise in the EW of both components of the doublet and a 
sharp decrease in the doublet ratio to 1.3 : 1. While the EW of the Fe II line also decreases 
throughout the first three spectra, it seems to completely disappear in the fourth spectrum. 
The two small features in this spectrum at the location of the Fe II feature are consistent 
with noise (see Figure [2]). The Fe II absorption line gradually returns throughout the last 
three spectra, but it never regains its original strength. The only other spectral feature to 
show possible variability is the Si IV A1393.76 absorption line at the presumed redshift of the 
host galaxy. The variability in this line does not follow the behavior of either the Mg II or 
Fe II systems at z = 1.48. We are still investigating the nature and significance of variability 
in this line. 

Variation in the EW of an absorption line is a direct consequence of variation in the 
column density of the absorber as seen by the background source. The Fe II and Mg II sys- 
tems are not necessarily expected to be distributed homogeneously throughout an absorbing 
"cloud." As the beam size of the GRB increases due to its physical expansion, the beam 
could include various dense or dilute regions of the "cloud," leading to different behaviors in 
the measured EWs of the absorbers. Additionally, the Fe II and Mg II may not be co-spatial, 
so their behaviors may not be directly comparable as seen here. 

To estimate the statistical significance of these observed variations, we apply x 2 statistics 
to the "null hypothesis" that the Fe II and Mg II lines are not variable. Fitting the multiple 
measurements of the Fe II EW with a constant value gives \ 2 = 38.8 with six degrees of 
freedom. Such a value of \ 2 results in a probability of p — 0.0000794% (~ 5cr Gaussian 
significance) that the Fe II line is not varying. A similar analysis for the Mg II A2803 line 
results in x 2 — 14.0, which gives p = 2.93% (2.2a Gaussian significance) for lack of variation 
in the line EW. Finally, for the Mg II A2796 line, we find \ 2 — 7-4, resulting in p = 28.9% 
for the assumption that the EW does not vary. The combined data clearly show statistically 
significant evidence for variability. Note that while we only test here the significance of 
variations in the line EWs, it is also immediately apparent from the spectra that the line 
profiles are varying as well (see Figure [2]). Each of these lines is likely composed of several 
components that are blended in our spectra, and it is variations in individual components 
that are causing the observed variations in the total EWs and line profiles. 

To further investigate the variability of the lines, we use the spectral subtraction algo- 
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rithm developed by Hartman & Stanek (2007, in preparation), which uses the entire spec- 
trum in the wavelength range 6250A ~ 7200A to obtain the convolution kernel for matching 
individual observed spectra to the reference spectrum. We choose the first observed spec- 
trum as the reference, then calculate the difference between the matched spectra and the 
observed spectra for the rest of the observed spectra, divided by the variance spectra we 
can calculate the significance of the difference. The Fe II shows a 4.9a variation in the line 
region (6439 A ~ 6450A) and the Mg II A2796 line shows a 2.5a variation in the line region 
(6926 A ~ 6942A). 



4. Discussion 



Variable intervening Fe II and Mg II was predicted by iFrank et al. I (120061 ) as a direct 



consequence o f their geometri c solut ion to the Mg II absorption over-abundance problem 



discovered by iProchter et al. I (120061 ). Our result immediately provides constraints on the 
physical size of the Mg II absorbers that are several or ders of magnitude s tronger than those 

0.01 pc jFrank et al. Ibood ) vs. < 200-300 pc 



(Rauch et al. 


2002 


proposed by '. 


Dong 



cf. Petitjean et al. Il2000l ). Observations of lensed quasars such as those 



(120061 ) may provide additional constraints on the physical sizes of the 
absorbers. Significant spectral variability on the timescales of months to years induced 
by microlensing can be used to probe the 10 14 — 10 16 cm cloudlets along the sightlines of 
lensed quasars. Properties of the cloudlets (physical size, sub-structure, etc.) can be studied 
statistically if such spectral variations are seen. 

We can also investigate the average physical density of the absorbers. Assuming the 
"cloud" is optically thin, we can obtain a lower limit on the Mg II column density A^M g n = 
1.7 x 10 13 cm" 2 . To evaluate this result, we have used our third spectrum taken ~ 5 hours 
after the burst, where the Mg II line ratio is close to 2 : 1, indicating the lines are not 
strongly s aturated (see Table [IJ). The beam s ize of the GRB during that epoch can be 
estimated (jWaxman 1119971 : lLoeb fe Perna Ill998l ) via 



R s (t) = 4.1 x 10 



15 



E, 



52 



ni 



[l + z s )-^(t/hourfl % 



cm 



X 



where z s is the source redshift, £"52 is the "isotropic-equivalent" of the energy release in units 
of 10 52 er gs" 1 , and is the ambient gas density in units of 1 cm -3 . For GRB 060206, 
-E52 = 5.8 (jPalmer et al. 1120061 ). and assuming n\ = 1 cm -3 , then R s = 5.1 x 10 15 cm for our 



third spectrum. With the physical size of the individual Mg II "clumps" constrained to be 
comparable to the GRB beam size (i.e., ~ 10 16 cm), one can immediately derive the average 
physical density of Mg II atoms to be about 0.0015 cm -3 . This is the first time that such 
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a strong constraint has been obtained. Given its large Mg II equivalent width, the z=1.48 
absorber is likely a Lyman limit system with log N(HI) > 18, although there is a ~ 20 — 30% 
probability that it is a damped Lyman-alpha system (Rao et al. 2006). Then the physical 
density of the hydrogen of this cloud is > 200 cm -3 . The estimate of physical density of 
Mg II absorbers is very important for understanding the cloud origins and properties, and 
more generally the metal lines and metal enrichment (see discussion in Schaye et al 2007). 

Another consequence of this discovery relates to GRBs themselves. In the case of 
GRB 060206, we see the Mg II EW initially decreasing, but then it shows a dramatic in- 
crease. This behavior is either a result of the GRB beam becoming much smaller in a short 
time (possibly due to an episode of energy injection) or part of the afterglow ring suddenly 
brightening. Our discovery points the way to a new method for studying GRB phenomena 
by monitoring afterglows with high-resolution optical spectroscopy. 

Additional implications relate to chemical abundances for intervening systems as mea- 
sured using absorption lines. In our first spectrum, the ratio of EWs for Fe II compared 
to Mg II is about 0.7. In our fourth spectrum, where the Fe II is only marginally present, 
and Mg II is very strong, that ratio is less than 0.1. Estimates of relative abundances for 
these two elements would be dramatically different when analyzing these two low-resolution 
spectra. If the system is at the lower end of this column density rangelog N(HI) > 18, some 
of the difference in Mg II/Fe II ratio could be the result of spatially- varying ionization lev- 
els. However, it may be that there are intrinsic variations in the relative metal abundances 
(although such variations are not seen in damped Lyman-alpha systems: Prochaska 2003, 
Rodriguez et al 2006). 

The observations presented in this paper open up the exciting possibility of study- 
ing time-variability of intervening absorption lines in the spectra of cosmological objects. 
Strongly and rapidly variable objects, such as GRBs and blazars, are obvious targets for fur- 
ther studies, but quasars (including both lensed and non-lensed) should also be investigated. 
Variable narrow absorption lines have been observed in spectra of quasars (e.g., Hamann et 
al. 1997), but they have been interpreted as originating in gas physically associated with 
the quasar. In fact, variability of narrow absorption lines is now used as one of the criteria 
to select intrinsic absorbers associated with quasars. Results presented in this paper suggest 
that using this criterion should be re-examined — quasars are optically variable sources, so 
size of their continuum emission region is bound to vary, just like it does for GRBs. As 
already discussed by Frank et al. (2006), quasars with smaller continuum emitting regions 
(i.e., intrinsically fainter ones), should have higher incidence of the strong Mg II absorbers. 

Note that throughout this paper we considered the z = 1.48 variable absorption system 
not to be physically associated with the z = 4.05 GRB, i.e., to be truly intervening due to 
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cosmological expansion. For these two objects to be physically associated would require Fe II 
and Mg II absorbing gas to have ejection velocity v ~ 0.6 c, while maintaining relatively low 
temperature in order not to be completely ionized. We consider such scenario most unlikely. 
Detection of a galaxy associated with this particular z = 1.48 variable absorption system 
would settle this issue beyond any reasonable doubt. As discussed in Prochter et al. (2006), 
such foreground galaxies were already identified for a number of strong, intervening Mg II 
absorbers seen in GRB spectra. 

We thank J. Hartman, P. Martini, S. Gaudi, S. Mathur, R. Pogge, S. Frank, X. Dai, 
J. Prieto, D. An, M. Elvis, R. Kirshner, and R. Narayan for helpful comments and discussions. 
M. C. Bentz is supported by a National Science Foundation Graduate Fellowship. 
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Table 1. EVOLUTION OF EQUIVALENT WIDTH OF THE z = 1.48 ABSORPTION 

LINES 



tobs 


Age b 


EW C Mg II 2796 


EWMg II 2803 


EWFe II 2600 


2.86384 


4.13 


2.53 ± 0.23 


2.05 ± 0.25 


1.72 ± 0.25 


2.88946 


4.50 


2.14 ± 0.25 


1.50 ± 0.27 


1.30 ± 0.15 


2.91110 


5.02 


2.06 ± 0.30 


0.92 ± 0.33 


0.60 ± 0.20 


2.93323 


5.55 


3.14 ± 0.44 


2.47 ± 0.41 


0.28 ± 0.21 


2.95857 


6.15 


2.24 ± 0.29 


1.93 ± 0.34 


0.36 ± 0.17 


2.99745 


7.09 


1.90 ± 0.44 


2.34 ± 0.58 


0.96 ± 0.29 


3.01990 


7.63 


1.96 ± 0.35 


1.41 ± 0.42 


0.67 ± 0.40 



a Helio centric Julian Date at middle of exposure minus 
2,453,770. 

b Age of the middle of exposure in h ours from the detect ion of 
the burst at 2006 Feb 06 04:46:53 UT jMorris et al. Il2006h . 

c Observer-frame equivalent width. 
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Fig. 1. — Spectra of the OT associated with GRB 060206. Prominent absorption features 
associated with three redshift systems are indicated: z = 1.48 (black long dashed vertical 
lines); z = 4.05 (red solid vertical lines); z = 3.54 (blue short dashed vertical lines). 
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Fig. 2. — Spectra of the GRB 060206 taken at 4.13 hours, 5.02 hours and 5.55 hours after 
the burst trigger in the z = 1.48 Fe II and Mg II region. Each spectrum has been offset by 
a constant flux value for clarity. 
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4 5 6 7 

Time after the burst [hours] 

Fig. 3. — Time evolution of the observed EW of the absorption lines of Mg II A2796 (filled 
circles) Mg II A2803 (open circles), and Fe II A2600 (filled squares). The EW measurements 
for the different lines have been slightly offset in time for clarity. 



